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■ We study the phase-ordering kinetics following a quench to a final temperature Tf of the one- 
dimensional p-state clock model. We show the existence of a critical value p c = 4, where the 

f») | properties of the dynamics change. At Tf = 0, for p < p c the dynamics is analogous to that of the 

. kinetic Ising model, characterized by Brownian motion and annihilation of interfaces. Dynamical 

' scaling is obeyed with the same dynamical exponents and scaling functions of the Ising model. For 

p > pc, instead, the dynamics is dominated by a texture mechanism analogous to the one-dimensional 
XY model, and dynamical scaling is violated. During the phase-ordering process at Tf > 0, before 
equilibration occurs, a cross-over between an early XY-like regime and a late Ising-like dynamics is 
observed for p > p c . 
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I. INTRODUCTION 



> 

CO 

After quenching a ferromagnetic system to a low temperature phase, relaxation towards the new equilibrium state is 
realized by a progressive phase-ordering pj . The specific mechanisms involved in the coarsening phenomenon depend 
3 ■ on the presence and on the nature of topological defects seeded by the disordered initial configuration which, in turn, 
are determined by the space dimensionality d and the number of components N of the order parameter. For N < d 
defects are spatially extended; in this case coarsening is driven by reducing the typical curvature of the defect core, 
removal of sharp features and shrinking of domain bubbles or vortex loops. Systems with N = d are characterized 
by the presence of stable localized topological defects and the ordering process occurs by mutual defect-antidefect 
annihilation. This is the case of the Ising chain quenched to a final temperature Tf = 0, where up and down domains 
are separated by point-like interfaces performing Brownian walks. When N = d + 1, such as in the one dimensional 
^ XY model, the kinetics is characterized by textures, spatially extended defects without a core, along which the order 
^C) ] parameter rotates by 2n. Growth of the typical size of textures is a relevant mechanism at work in these systems. 
Finally, for N > d + 1 topological defects are unstable and the dynamics is solely driven by the reduction of the excess 
' energy related to the smooth rotations of the order parameter. 

In any case, the development of order is associated to the growth of one or more characteristic lengths, with laws 
that, besides the specific mechanisms discussed above, depend on the conservation laws of the dynamics. 

Generally, the late stage is characterized by dynamical scaling. This implies that a single characteristic length 
L(t) can be associated to the development of order in such a way that configurations of the system are statistically 
independent of time when lengths are measured in units of L(t). The characteristic length usually has a power law 
■ growth L(t) oc t x l z . In systems with a non-conserved order parameter one generally finds z = 2. In particular, this 
value is provided by the exact solution of the kinetic Ising chain quenched to zero temperature. 

However, there are cases where dynamical scaling is violated, notably the XY model in d = 1,2. In d = 1 this 
is related Q to the presence of two lengths L w (t) and L c (t), associated to the texture length and to the texture- 
antitexture distance, growing with different exponents z — 4 and z — 2 respectively. 

In this Article, we investigate the interplay between two coarsening mechanisms, point-like defect annihilation and 
texture growth, in the phase-ordering kinetics of the one dimensional p-state clock model. This spin system reduces 
to the Ising model for p = 2 and to the XY model for p = oo. We study how the model with generic p interpolates 
between these limiting cases which, as discussed above, behave in a radically different way. In doing that, we uncover 
the existence of a critical value p c — 4, where the properties of the dynamics change abruptly. For p <p c the dynamics 
at Tf = is characterized by Brownian motion and annihilation of interfaces between domains, as in the Ising model. 
One has dynamical scaling with the same dynamical exponents and, interestingly, the same scaling functions of the 
Ising model. For p > p c , instead, the dynamics is dominated by a texture mechanism analogous to the case with 
p = oo, and dynamical scaling is violated. 

In d = 1 there is no possibility of ergodicity breaking except at T = 0. At any finite temperature the equilibrium 
state is disordered with a vanishing magnetization and a coherence length £(T) that diverges in the T — > limit. If the 
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system is quenched to a sufficiently low temperature one has a coarsening phenomenon in a pre-asymptotic transient 
until the growing length associated to the development of order becomes comparable with £(7/). Since £(T/) diverges 
as Tf — > the phase-ordering stage can be rather long. In this regime we show that activated processes restore, after 
a characteristic time T£ ross (Tf), the Ising behavior also in the cases with p > p c . 

This paper is organized as follows: In Sec. [H] we introduce the model and define the observable quantities that 
will be considered. In Sec. II I II we present the outcome of numerical simulations of the model with different p. In 
particular, quenches to Tf = or to Tf > will be discussed in Sees. IIII Al and ITlI Bl respectively. Here we enlighten 
the crossover between the Ising and the XY universality class and provide an argument explaining its microscopic 
origin. A summary and the conclusions are contained in Sec. IIVI 



II. MODEL AND OBSERVABLES 



The p-state clock model in one dimension is defined by the Hamiltonian 



H[a] = -J^<? ( ■ ffi+i = -J^2cos(6i - i+1 ), (1) 

i=l i=l 

where is a two-components unit vector spin pointing along one of the directions 

9i = —n h (2) 
P 

with rii G {1, 2, ...,p}, i = l,...,J\f are the sites of the lattice and we assume periodic boundary conditions 0^+1 = 
This spin system is equivalent to the Ising model for p = 2 and to the XY model for p — > oo. In d = 1 the system is 
ergodic except at T = 0. At any finite temperature the equilibrium state is disordered with a vanishing magnetization 
and a coherence length £(T) that diverges in the T — > limit. 

We consider a system initially prepared in an high temperature uncorrelated state and then quenched, at time t = 0, 
to a lower final temperature Tf. The dynamics is characterized by the ordering of the system over a characteristic 
length growing in time until, at time Tp q (Tf) it becomes comparable to £(?/). At this point the final equilibrium 
state at Tf is entered. Quenching to Tf — 0, since £(0) = oo, one has Tp q (Tf) = oo; therefore an infinite system never 
reaches equilibrium and the phase-ordering kinetics continues indefinitely. If the system is quenched to a sufficiently 
low temperature, since is very large, the same behavior, as for Tf = 0, can be observed over the time window 

t<T p q{ ?f)- 

The power law growth of the characteristic size of ordered regions depends on the specific mechanisms at work 
in the kinetic process. In the Id Ising model with non conserved order parameter, i.e. single spin flip dynamics, 
ordering is determined by the Brownian motion of the interfaces between up and down domains, which annihilate 
upon meeting. This leads to 

L(t)~i-, (3) 

with z = 2. The same value is also expected Q for p < 4. 

The situation is different in the XY model in d = 1. Here the order parameter is a vector which can gradually 
rotate with a low energy cost. A smooth 2tt rotation of the phase 8 is called a texture when the rotation is clockwise, 
or antitexture when it is counterclockwise. The length over which this phase winding occurs will be denoted by L w {t). 
After a quench from a disordered state textures and antitextures are formed with equal probability. Then, there are 
points where the rotation of 9 changes direction and the phase decohere. We denote with L c (t) the characteristic 
length over which the phase remains coherent. It was shown [3j that L w (t) and L c (t) grow with a power law @ but 
with different exponents. Specifically one has z = 4 for L w (t) and z = 2 for L c (t). The existence of these two lengths 
is at the heart of the scaling violations of the XY model. 

Characteristic lengths can be estimated from the knowledge of the two-points equal time correlation function 

G(r,t) = (ai(t)-a i+r (t)), (4) 

where (. . .) means an ensemble average, namely taken over different initial conditions and thermal histories. Due 
to space homogeneity, G(r,t) does not depend on i. If there is a single characteristic length in the system, one has 
dynamical scaling which implies 



G(r,t)=g(x), 



(5) 
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where x = r/L(t). In the Ising model one finds Q 

g(x) = erfc{x}. (6) 



with L(t) = \j2t. For small x one has the Porod linear behavior 1 — g(x) ~ x, which is expected for systems with 
sharp interfaces 1]. From Eq. JSJ one can extract a quantity Lq(J,) proportional to L(t) from the condition 

G(L G (t),t) = i (7) 

namely as the half-height width of G(r,t). In the XY model, G(r,t) still obeys Eq. JSJ, with x = r/L w (t), L w [t) — 
2 3 / 4 (^) 4 / 4 , and Q 



5 (a;)=exp|- — |. (8) 

Here £j is the correlation length of the initial condition which, for a quench from a disordered state, is of the order of 
the lattice spacing. The Porod law is not obeyed, since instead of sharp interfaces one has smooth textures. Note that 
G(r, t) has a scaling form, although dynamical scaling is violated. Scaling violations can be evidenced by considering 
different quantities as, for instance, the autocorrelation function 

C(t,s) = (a i (t)-a i (s)). (9) 

In the Ising model this quantity can be cast in scaling form 

C(t, 8 ) = h(y), (10) 

where y — t/s and 



2 / 2 

Kv) = -arcsinW— — . (11) 
7r v 1 + y 

In the XY model, instead, one finds the stretched exponential behavior 

C(t, s) = cxp [-^s 1 [2(2/ + 1)* - (22/)* - V2] J , (12) 
This expression cannot be cast in a scaling form, as for the Ising model, revealing the absence of dynamical scaling. 



III. NUMERICAL RESULTS 



In the following we will present the numerical results. Setting J = 1, for each case considered we simulated a string 
of 10 4 spins with periodic boundary conditions and different values of p ranging from p = 2, corresponding to the Ising 
model, to p — 00, corresponding to the XY model. We consider a single spin flip dynamics regulated by transition 
rates 

w{[a] - [a']} = = I exp(^) 

T p exp(^) +exp(^^) 

Here [a] and [a'] are the spin configurations before and after the move, differing at most by the value of the spin 
on a randomly chosen site, AE = H[a'] — H[a], and we have set the Boltzmann constant to unity. The transition 
rates (|13f) are a generalization of Glauber transition rates to the p-state spins of the clock model. They reduce to 
the usual Glauber transition rates w{[a] -»■ [a 1 ]} = (1/2)[1 + tanh(-AE/2T)] for p = 2. The factor 2/p in Eq. JT3J 
ensures that all spin values have the same probability 1/p when AE = 0. 

An average over 10 4 realizations is made for each simulation. The statistical errors in the data reported in the 
figures are always smaller than the dimension of the symbols or the thickness of the lines. 
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A. Quenches to Tf = 0. 

Let us start with quenches to Tf — 0, by illustrating the behavior of the characteristic length La(t) defined in 
Eq. Q. In Fig. ^this quantity is plotted against t 1 / 2 (left panel) or against t 1 / 4 (right panel), for several values of p 
ranging from p — 2 to p — oo. 

This figure shows that Lc(t) has an asymptotic power law growth, as in Eq. ©, for every value of p. However, 
the dynamic exponent z radically changes going from p < p c , where one has values very well consistent with z = 2 
(best fits yield 1/z — 0.49 ± 0.01 for p = 2,3,4), to p > p c where z — 4 is found with good accuracy (we find 
l/z = 0.27 ± 0.0f,0.27±0.0f,0.25± 0.01 for p = 5,6, 10). We recall that these are the values found in the Ising model 
and in the XY model. The behavior of Lq{£), then, indicates a crossover from Ising to XY behavior upon crossing 
p c = 4. We will see in the following that this is confirmed by the analysis of other dynamical quantities. Before doing 
that, however, let us discuss which is the microscopic mechanism at the basis of this crossover. 

For finite values of 2 < p < oo we generalize the definition of a texture as a region of the lattice of length L w (i) 
where p subsequent domains are found, each of average length Ld(t) ~ L w (t)/p, such that moving along the lattice 
the value of n follows the sequence n= p,p— 1, 1. This is schematically shown in Fig. [3 

To developed order two mechanisms are possible: A texture can grow by increasing the number of spins L c i(t) on 
every step. We anticipate that this process is found to be relevant for p > p c and leads to the power law behavior © 
of L w (t) with z — 4, as in the XY model. This behavior competes with the tendency to build the largest possible 
domains, instead of textures. This amounts to replace a texture with a number Nd "C p of domains each characterized 
by a single value of n. However, for p > p c at T = 0, once textures are present, this process is not allowed. In fact, 
let us consider the situation of Fig. [2] and the possibility to form, in this region, a unique domain with, say, n = p 
(the dotted line in Fig. |2J) . There are several ways to do this. Suppose one starts by rotating the spins with n = 1 to 
n = 2, as shown by the thin arrow in Fig. [21 After the move the energy would change by an amount 

AE P = J[2cos(2tt/p) - cos(4tt/p) - 1]. (14) 

This function is plotted in Fig. |3J Interestingly one has AE p < or AE p > for p < p c or p > p c , respectively. 
At Tf = moves with AE p > are forbidden. Therefore, for p > p c there is no possibility to destroy the textures 
and form domains. Other possible moves, as, for instance, a rotation from n = 1 to n = 3, correspond to a larger 
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FIG. 2: (Color online) Schematic representation of a generalized texture for p — 5. 
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FIG. 4: (Color online) The average length of textures is plotted against t 1 



activation energy and are forbidden as well. Therefore, for p > p c textures and antitextures are stable against domain 
formation and the only ordering mechanism left is their growth and annihilation, much in the same way as in the XY 
model, leading to z = 4. Conversely, for p < p c textures are removed and domains are created whose competition 
leads to the Ising like behavior z = 2. As already discussed, in the XY model the exponent z = 4 is associated 
to the growth of the size of single textures. In order to check if the same mechanism is at work also in the clock 
model, in the numerical simulation we have identified the textures present in the system at each time and we have 
computed their average size L w (t). The results are shown in Fig. 0] for different values of p > p c , showing that, 
actually, the size of textures grows as a power law L w ~ t 1 ^ with z quite compatible with z — 4 (best fits yield 
l/z = 0.29 ± 0.02, 0.29 ± 0.02, 0.28 ± 0.02, 0.23 ± 0.02 for p = 5, 6, 10, 25, respectively). This confirms that the exponent 
z = 4 of the algebraic growth of La(t) is determined by the texture mechanism, as in the XY model. 

The previous results for Lc(t) indicate the presence of a crossover at p = p c from the Ising to the XY non- 
equilibrium universality class. In order to substantiate this conjecture we have computed other dynamical quantities. 
The equal-time correlation function is plotted in Figs. 1516171 against x — r/Lc(t). In Fig.[5]the cases with p = 2,3,4 
are considered. According to Eq. (JjJJ for p = 2 one should find collapse of the curves with different s on a single 
mastercurve g (x) given by Eq. 10 . This is indeed observed in Fig. El According to our hypothesis the same behavior 
should be observed also for p — 3,4, as can be verified in the figure. Moreover, one also finds that the mastercurves 
g(x) are numerically indistinguishable for different p, and they all coincide with that of Eq. JSJ. This result is trivial 
for p = 4, since in this case the clock model can be mapped exactly on two non-interacting Ising models. The same 
property could be expected also for p — 3. In fact, by considering G(r,i), it easy (see Appendix) to check that 

G(r,t) = ^G P (r,t)~±, (15) 

where Gp(r,t) is the single phase equal time correlation function of the 3-state Potts model. This quantity was 
computed in J5|, where it was found 

G P {r,t) = ^G I [r,t) + 1 -, (16) 
where Gi (r, t) is the equal time correlation function of the Ising model. Plugging Eq. (|16|) into Eq. I|15|l one finds 
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FIG. 5: (Color online) The correlation function G(r,t) is plotted against x = r / La(t) for p — 2, 3, 4 at t = 1800. The dashed 



G(r, t) = Gi(r, t). The same argument shows also the identity between the two time correlation functions of the clock 
model with p = 3 and the Ising model, strongly suggesting the complete equivalence between these models. 

Let us emphasize that this result indicates a stronger similarity among the cases p = 2,3,4 than a unique non- 
equilibrium universality class would imply, since not only the exponents are equal but the whole functional form of the 
scaling function. This results are in contrast with those of ref. la where an approximate theory was used to show the 
dependence of g(x) on p. However, the approximation used in jy| is expected to improve increasing the dimensionality 
d. 

The cases with p > p c are shown in Fig. 16171 As discussed in Section [H] G(r,t) obeys the scaling form JSJ) also in 
the XY model, although dynamical scaling is violated. According to our conjecture, for p > p c we expect the same 
behavior. In Fig. [51 it is shown that, indeed, the curves at different times collapse when plotted against x = r/La(t). 
However, differently from the cases p < p c , the masterfunction g(x) depends on p and converges to the form (JSJ of 
the XY model for p — ► oo, as shown in Fig. [7| 

Let us turn to consider the autocorrelation function, that is plotted in Figs . 181 lUl against y — t/s. In Fig.|S|the cases 
with p = 2, 3, 4 are considered. Here the situation is analogous to that of G(r, t). For p = 2 one should find collapse of 
the curves with different s on a mastercurve h(y), Eq. IjlOl) . This is indeed observed in Fig. EI The same behavior is 
observed also for p — 3,4. Again, as for G(r, t), we find that the mastercurves h(y) are numerically indistinguishable 
for different p, and they all coincide with that of Eq. (jl 1|> . In order to check if this property is completely general, 
namely if every observable is characterized by the same exponents and scaling functions for p = 2,3,4, besides the 
correlation functions we have also computed the integrated autoresponse function 



line is the analytic expression JSJ. 




(17) 



Here 




(18) 



a = 1,2 being the generic vector components, is the linear autoresponse function associated to the perturbation 
caused by an impulsive magnetic field hi switched on at time t' < t. In the Ising model 0, in the T — > limit one 
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FIG. 6: (Color online) Data collapse of the correlation function G(r,t) plotted against x — r/La(t) for p — 5 (upper panel) 
and p = 6 (lower panel), at different times (t = 190, 245, 315, 405, 520, 665, 855, 1100, 1400, 1800). 
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FIG. 7: (Color online) The correlation function G(r,t) is plotted against x = r/La(t) for different values of p, at t = 1800. 
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FIG. 8: (Color online) The autocorrelation function is plotted against y for p = 2, 3, 4. The dashed line is the analytic 
expression 1111 . 



finds 

X (t,s) = f(t/s), (19) 

with 

f(y) = -L(l-?- arcsin y/^j . (20) 

Here we measure the response function using the efficient method derived in [8j without applying the perturbation. 

The behavior of x(t> s ) is shown in Fig. [21 for the cases p = 2,3, 4. One finds collapse of the curves with different s 
on a mastercurve f(y), as in Eq. (|19|l for p = 2. Also in this case mastercurves f(y) for different p are numerically 
indistinguishable. In conclusion, then, our data for G(r,t), C(t,s) and x(^: s ) confirm that the cases with p = 2, 3,4 
share the same exponents and scaling functions. Notice that having the same scaling function both for C(t,s) and 
x(t, s), the cases with p < p c have also the same parametric plot of x(t, s) versus C(t, s) 0. 

The situation is radically different for p > p c . We expect here to see a texture-dominated XY-like dynamics, with 
violations of dynamical scaling that can be detected from C(t, s). In fact, this is what one observes in Fig. where 
the autocorrelation function is plotted against y. For each value of p, curves with different values of s do not collapse. 
The whole behavior is qualitatively similar to that of the XY model described by Eq. i|12H , which predicts the lowering 
of the curves for fixed y as s increases. Quantitatively, as already observed regarding G{r,t), the analytic form of the 
curves depends on p and is different from that of the XY model, namely Eq. I|12|) . As shown in Fig. ^2 Eq. I|12l) is 
gradually approached increasing p. 



B. Quenches to T f > 0. 

When quenches to finite temperatures are considered, as already discussed in Sec. [HI one has a finite equilibration 
time Tp q {Tf). In the following we will always discuss the ordering kinetics preceding the equilibration time, namely 
for t <C Tp q (Tf). 
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FIG. 9: (Color online) s ) i s plotted against y for p — 2,3, 4. The dashed line is the analytic expression 12011 . 



According to our hypothesis, the XY-like behavior observed for p > p c is due to the impossibility to eliminate 
textures and form domains, because this would require activated processes with AE p > given by Eq. (|14|l . Quenching 
to a finite temperature those processes are no longer forbidden and we expect textures to start being removed after 
a characteristic time Tp ross (Tf). In order to estimate the crossover time let us consider again the situation of Fig. [3] 
The activated process described by the thin arrow, where the spins with n = 1 are rotated to n = 2, is a first action 
towards the removal of the texture, but the texture is not disappeared yet. The second action is the rotation of spins 
from n — 2 to n — 3, indicated by a bold arrow in the figure |9(. This requires an energy 

AE p 2) = J[cos(2tt/p) + cos(4tt/p) - cos(6tt/p) - 1]. (21) 

Then a third action is required, where spins with n = 3 are rotated to n = 4, and so on, until, after p — 1 steps all 
the spins in the region considered have n = p. It is easy to generalize Eqs. <|14I21H to the generic m-th action: 

AE^ = J[cos(2tt/p) + cos(2m7r/p) - cos(2(m + l)ir/p) - 1]. (22) 

(2) 

Let us consider AEp . This quantity is positive for p > 6. For p = 5,6, therefore, the second action is not an 
activated process, while it is activated for p > 6. In general, from Eq. H22[) one has AE { p m) > for p > 2 + 2m. The 
accomplishment of an action requires a time |10| 

t p m \T f ) ~ [wpiAEp/TfT 1 = I {l + cxp[A4™>/7>]} , (23) 

w p being the transition rates defined in Eq. (|13|) . The crossover time, namely the characteristic time after which 
textures are removed, is given by the sum of the times required for all the p — 1 actions. It can be evaluated as 

p-i 

T cross [Ts)= Y j tp m \Tf). (24) 
m=l 

In the limit Tf — > the sum is dominated by the process with the largest activation energy 

T cross [Tf „ Q) = SuP{m=1 , v _ 1} t p m \T f ). (25) 

The Sup in this equation is obtained for m = ni* given by 

1 for p < 10 



[£z2] for p > 10 



(26) 
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FIG. 10: (Color online) The autocorrelation function is plotted against y for p — 5 (upper panel) and p = 6 lower panel). 



where [x] is the integer part of x. Then, in the low-T limit one has 

Tp ross (Tf ~ 0) = t™* (Tf). (27) 

In conclusion, for p < p c no activated processes are required and the system immediately enters the Ising-like phase 
ordering behavior. For p > p c , instead, the dynamics is initially of the XY type until, at t ~ T p TOSS (Tf) there is a 
crossover to the Ising-like non-equilibrium behavior. 

The crossover can be appreciated in Figs. II 21131 The former shows the behavior of Lc(t) for p = 6 and different 
values of Tf. Here one observes initially the same behavior as for Tj = 0, namely Lc{t) oc W 4 , i.e. a straight line in the 
plot of Lait) against i 1 / 4 (right panel). For larger times there is a crossover to the Ising behavior Lq{£) ~ t 1 / 2 , namely 
a straight line in the plot of Lc{i) versus t x l 2 (left panel). Although the crossover is a quite smooth phenomenon, as 
can be seen in Fig.flEZl T^ ross {Tf) given by Eq. J23J, represented by thick segments across the lines, turns out to be 
of the correct order of magnitude for all the temperatures considered. 

In Fig. ^5] we plot Lc(t) for Tf —0.1 and different values of p. One observes the same pattern of behavior of Fig. 1121 
with a crossover from a power law growth with z = 4 to one with z = 2. The crossover time (|24|l grows with p, as 
expected. 



IV. CONCLUSIONS 



In this paper we have studied the phase-ordering kinetics of the one dimensional p-state clock model. We have 
shown the existence of a critical value p c — 4 separating two radically different dynamical behaviors. For p < p c 
the dynamics is in all respects analogous to that of the Ising model with p = 2. Phase-ordering proceeds by means 
of formation and subsequent growth of domains through interface diffusion and annihilation. This similarity goes 
beyond the qualitative level: we find the same exponent and scaling functions for every p < p c and for all the one-time 
or two-time quantities considered. This reflects a deeper similarity than what a unique universality class, involving 
only the value of the exponents, would imply. For p > p c the dynamics changes dramatically, due to the relevant 
role played by textures. While for p < p c textures are quickly removed by means of non-activated processes, for 
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FIG. 11: (Color online) The autocorrelation function is plotted against y for s — 100 and several values of p. 




FIG. 12: (Color online) Le(i) is plotted against t 1 ^ 4 (left panel) or versus t 1 ^ 2 (right panel), for a quench of a system with 
p = 6 and different values of Tf. Vertical segments on the curves for different T> represent Tp roaa (T f ) obtained from Eq. 
as discussed in the text. For the smallest temperature Tp roas (Tf) is outside the range of times of the figure. 



p > p c their removal can only be realized through activated processes. For quenches to Tf = 0, activated process are 
forbidden, and, therefore, textures remain in the system up to the longest times. Their peculiar growth mechanisms 
characterize the dynamics, similarly to what happens in the one-dimensional XY model, with the notable feature of 
violation of dynamical scaling and the anomalous growth with z — 4 of the winding length L w (t). For quenches to 
finite Tf, textures survives up to a characteristic time Tp TOSS (Tf) which can be rather long for small temperatures 
or large p. A crossover phenomenon is then observed from an initial dynamics of the XY type, to a later Ising-like 
behavior. 

Our results are at odd with what is found in Ref. ja] where an approximate analytical solution of the clock model 
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FIG. 13: (Color online) -Lg(£) is plotted against t 1 ^ 4 (left panel) or versus t 1 ^ 2 (right panel), for a quench at Tf = 0.1 and 
several values of p > p c . Vertical segments on the curves for different p represent Tp ro3S (Tf). For p > 25 Tp roas (Tf) is outside 
the range of times of the figure. 



in arbitrary dimension is obtained, finding an analogous scaling behavior for all p < oo but with p-dependent scaling 
functions. In the present one-dimensional case, instead, the situation is the opposite. There is not an analogous scaling 
behavior for all values of p, but a qualitative difference occurs crossing p c . In addition, when scaling holds, namely 
for p < p c , the scaling functions do not depend on p. We believe, however, the behavior of the system considered in 
this paper, to be peculiar. Actually, the different dynamics observed crossing p c is determined by the simultaneous 
presence of interfaces and textures. On the basis of the discussion of Sec. U] we expect a similar situation to be only 
realized in iV-component vectorial models with discrete states and N = d + 1, where extended defects without a 
core may exist. For instance, it would be very interesting to study if a similar pattern is observed in d = 2 for a 
generalization of the clock model where a three component order parameter is only allowed to point on a finite number 
p of directions. In addition, we expect the remarkable feature of unique scaling functions for different values of p to 
be peculiar to the one-dimensional case. Considering the function G(r, i), for instance, the scaling function describes 
the spatial distribution of domains and it is quite evident that in d > 1 this depends on p. Taking the case d = 2, for 
simplicity, one has the usual bicontinuous domain structure of domains and interfaces for p — 2, while for p > 2 there 
is a different pattern with interfaces and vortices |ll| . However, in the one dimensional case interfaces are point-like 
objects for all values of p and one does not expect relevant differences in their spatial distribution when p is changed. 

Finally, it would be very interesting to study if a similar pattern is observed in the one-dimensional clock model 
with a conserved order parameter. Concerning the value of the growth exponent z, which in the non-conserved case 
considered here effectively discriminate the Ising dynamics with z = 2 from the XY behavior with z — 4, in the 
conserved case one should observe a crossover from z = 3 to z = 6 0, • 
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APPENDIX 



For the 3-states clock model the correlation between two spins at a certain time t can be written as 
G{r,t) = (o-itrj) = coa M n ) - Oj{n')] p i( n ^) p iA n ^ I 



(28) 
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where r is the distance between i and j. n, Bi (and their relation) are defined in Eq. Pj(n, t) is the probability to 
find the spin on site i in the state n at time t, and Pi^(n,t \ n',t) is the conditional probability to find the state n' 
on site j provided that the state n is found in i. Isolating the diagonal terms one has 

G(r,t) = Pi(n,t)P id (n,t \ n,t) - ^ £ i>-(n,t) £ P id (n,t \ n',t), (29) 

n— 1,3 n— 1,3 n'^n 

where we have used the value cos(^ — 0j) = —1/2 when 6>,; ^ 0j. Since ^ n i^ n Pi,j(n,t \ n',t) = 1 — Pi : j(n,t \ n 7 t) 
one has 

G(r, i) = - i P ^ n > *) + 5 P *( n ' *) Pi J ( n » * I n > *) = _ \ + | S Pi ( n > (», * I ». *) (30) 

n— 1,3 n— 1,3 n— 1,3 

Let us turn now to the Potts model where a generic spin on site i can be found in the states labeled with m,i = 1, 2, 3. 
Following Ref. 0, we define an auxiliary field <f>i(n) such that (f>i(n) — 1 if m, = n, where n is a reference state, and 
4>i(n) = otherwise. The correlation of the auxiliary field is the single phase correlation function of the Potts model 
and can be written as 

G n {r,t) = (^(n)^-(n)) = Pi(n,t)P^(n,t \ n,t), (31) 

where the probabilities are defined analogously to the those of the clock model introduced above. Recognizing G n (r, t) 
in the last term of the right hand side of Eq. (|30|l one arrives at 

G{r,t) = -\ + \ J2 G ^)- ( 32 ) 

n=l,3 

Because of the rotational symmetry one has Gp{r, t) = G n (r, t) for all values of n and then one recovers Eq. i|T5)l. 
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